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ABSTRACT 

The  development  of  highly  sophisticated  and  complex  ejection  seats  to 
provide  safe  egress  from  aircraft  during  low  altitude,  high  speed  and  adverse 
attitude  ejections  addresses  the  need  for  a  reliable  computer  simulation 
capability  to  determine  the  trajectory  of  the  seat  and  occupant  under  these 
conditions.  The  occupant  response  to  increasingly  severe  accelerations 
profiles  must,  therefore,  be  addressed  if  the  simulation  model  is  to  be  able 
to  evaluate  the  behavior  of  the  seat/occupant  system.  Traditionally,  the  scat 
and  occupant  have  been  treated  in  combination  as  a  rigid  body  for  purposes  of 
trajectory  analysis.  However,  experience  with  ejection  seat  tests  has 
demonstrated  a  considerable  amount  of  relative  motion  between  the  occupant 
(typically  an  anthropomorphic  dummy)  and  the  seat  itself.  This  paper  presents 
a  lumped  mass,  spring  damper  mathematical  model  to  simulate  the  motion  of  the 
occupant's  C.G.  with  respect  to  the  seat  under  a  dynamic  ejection  environment. 
The  analysis  of  anthropomorphic  dummy  results  and  of  computer  generated 
biodynamic  simulation  data  used  in  the  evaluation  of  th«  model  will  be 
discussed. 

This  work  was  sponsored  by  the  Naval  Air  Systems  Command  (AIR-340B)  under 
Airtask  No.  F4 1400000. 
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INTRODUCTION 

The  relationship  between  the  seat/occupant  center  of  mass  and  the  rocket 
thrust  line  has  basn  of  great  concern  in  the  ejection  seat  community  ever 
since  the  introduction  ot  the  booster  rocker.  This  relationship  is  critical 
in  determining  the  motion  of  the  seat/occupant  during  rocket  thrusting,  and 
ultimately  s  ongly  affects  the  trajectory  achieved  during  ejection.  The 
determination  of  the  seat /occupant's  C.G.  under  static  or  1  G  conditions  is  s 
relatively  simple  procedure  and  the  variation  in  C.G.  location  for  different 
anthropomorphic  and  human  occupants  is  well  documented  (1,2).  However , 
experience  with  ejection  seat  tests  with  both  anthropomorphic  dummy  and  human 
occupants  has  demonstrated  a  considerable  amount  of  motion  between  the 
occupant  and  the  seat  itself,  and  attempts  have  been  made  to  quantify  and 
describe  the  resulting  C.G.  shift  (2,3). 

Computer  jimulations  of  ejection  seats  can  b*.  used  ,  t :i actively  to 
evaluate  the  trajectories  of  the  seat  and  occupant  during  low  altitude,  high 
speed  and  adverse  attitude  ejections  for  which  actual  track  tests  are  not 
feasible.  The  motion  of  the  occupant's  center  of  mass  and  the  resulting 
seat/occupant  C.G.  shift  caused  by  these  increasingly  severe  ejection 
conditions  must  be  taken  into  account  in  order  for  the  simulation  model  to 
effectively  evaluate  the  behavior  of  the  seat/occupant  system. 

This  paper  describes  the  investigation  undertaken  to  develop  and  evaluate 
a  mathematical  model  for  the  motion  of  the  center  of  mass  of  an  occupant  unde.r 
ejection  accelerations.  Since  the  vast  majority  of  ejection  seat  tests  use 
anthropomorphic  dummies,  the  investigation  first  addressed  dummy  C.G.  motion 
and  a  mathematical  Dynamic  C.G.  Model  was  developed  and  evaluated.  The  model 
was  then  exercised  using  human  biodynamic  simulation  data  in  order  to  assess 
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its  effectiveness  in  replicating  the  motion  of  the  center  of  mass  of  simulated 
human  occupants.  Finally,  the  model  was  incorporated  into  an  existing  6 
Deg.  >3-of-Freedom  trajectory  simulation  computer  program  which  is  presently 
being  employed  to  conduct  trajectory  analysis  for  the  Escapac  Replacement 
Program. 

METHODOLOGY 
Ejection  Tower  Tests 

The  standard  scientific  paradigm  was  used  in  this  investigation:  collect 
data,  analyze  data,  l'ora  hypothesis,  and  test  hypothesis.  No  tests  were 
actually  conducted,  but  rather  the  data  base  from  all  the  previous  tests 
conducted  on  the  Ejection  Tower  at  NADC  was  used.  The  only  criteria  used  to 
select  tests  from  this  data  base  for  further  investigation  were  that  the 
occupant  had  to  be  an  anthropomorphic  dummy,  and  both  seat  and  dummy  vertical 
accelerations  had  to  have  been  measured  and  the  oscillograph  records  were 
available.  These  criteria  resulted  in  52  candidate  tests  (see  Table  1).  For 
each  of  these  tests,  the  seat  and  dummy  vertical  accelerations  were  digitized 
and  stored  on  tape.  The  two  accelerations  were  then  subtracted  and  double 
integrated  (figure  1).  The  quantity  "d"  shown  in  figure  1  represents  the 
relative  displacement  between  the  seat  and  the  center  of  mass  of  the  dummy, 
and  is  in  fact  the  quantity  that  is  being  modeled.  There  are  some  serious 
problems  with  the  accuracy  of  this  calculated  relative  displacement.  First, 
the  seat  and  dummy  accelerations  are  very  similar,  so  that  the  subtraction 
removes  a  large  part  of  the  signal  and  leaves  the  noise  and  error  to 
constitute  a  more  major  portion  of  what  remains.  Secondly,  double  integration 
greatly  magnifies  small  errors,  and  is  particularly  sensitive  to  zero  biases 
such  as  can  occur  during  digitizing.  The  form  of  the  resulting  error  is  a 
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function  which  grows  with  time:  either  parabolic  or  cubic.  This  type  of  error 
is  evident  in  a  number  of  the  52  tests  but  it  is  not  easy  to  correct, 
especially  since  the  tests  are  not  recent  and  it  is  not  possible  to  go  back 
and  recalibrate  the  accelerometers,  for  instance.  Another  problem  with  the 
data  is  that  it  is  possible  that  the  dummy  did  not  remain  perfectly  vertical 
and  parallel  to  the  seat  back  throughout  the  test,  and  so  ;he  two 
accelerometers  being  compared  do  not  point  in  exactly  the  same  direction. 
Despite  all  these  inaccuracies  in  the  data,  information  can  still  be 
extracted.  For  any  given  test  one  can  never  be  sure  that  certain  features  of 
the  data  are  real  or  artifacts  of  the  errors,  but  when  all  52  tests  are 
examined  together  certain  facts  show  through. 

By  visual  inspection  of  the  C.G.  displacement  curves  it  was  possible  to 
extract  a  general  pattern  (figure  1)  and  certain  parameters  were  identified  as 
potentially  significant  and  these  were  collected  in  Table  1.  To  characterize 
this  shape,  the  times  and  displacements  at  the  first  and  second  turns  for  all 
52  tests  were  recorded,  along  with  other  parameters  which  conceivably  might 
determine  or  at  least  influence  the  shape  of  the  displacement  curve;  these 
parameters  were  entered  into  the  computer  for  plotting  and  statistical 
analysis.  Figures  2  and  3  show  scatter  diagrams  of  the  displacements  at  the 
first  and  second  turns  of  the  standard  displacement  curve.  It  must  be  noted 
that  the  displacement  at  the  first  turn  is  quite  consistent  across  tests,  but 
that  the  second  turn  is  much  more  variable.  This  agrees  with  the  preceding 
error  analysis,  which  pointed  out  that  the  error  should  increase  rapidly  with 
time.  It  should  be  noted  also  that  the  two  displacements  are  very  nearly  the 
same,  which  indicates  that  the  paradigm  should  be  trapezoidal.  Figure  4  shows 
the  times  of  the  two  turn3.  The  time  of  the  first  turn  is  rather  constant, 
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but  the  second  turn  seems  to  vary  with  the  test  conditions :  it  is  not  constant 
over  all  tests  and  its  variations  do  not  appear  random.  The  hypothesis  that 
the  displacements  at  the  two  turns  are  more  or  less  the  same  is  tested  in 
figure  5.  If  the  displacements  at  the  two  turns  were  in  fact  equal,  the 
points  would  fall  on  the  straight  line  in  the  figure.  The  data  is  not 
inconsistent  with  the  hypothesis  that  the  displacements  at  the  two  turns  are 
nearly  equal,  but  i*  certainly  doesn't  confirm  the  hypothesis  either.  The 
displacement  at  the  first  turn  is  then  plotted  in  terms  of  the  peak  seat 
acceleration  and  the  maximum  seat  velocity  in  an  attempt  to  discover 
controlling  factors  for  the  displacement  curve  (figures  6  and  7).  The  figures 
show  no  pattern  at  all,  indicating  that  the  magnitude  of  the  displacement  is 
indipendent  (within  the  small  range  examined)  of  the  strength  of  the  input 
acceleration.  Figures  8  and  9  test  the  hypothesis  that  the  second 
displacement  turn  is  caused  by  a  release  of  force.  The  release  of  force  was 
represented  by  the  time  of  thrust  end  (which  was  sometimes  difficult  to 
meacure)  (figure  8),  and  the  time  of  maximum  seat  velocity  (i.e.,  the  time 
when  the  acceleration  crosses  zero  and  goes  negative)  (figure  9).  In  both 
cases,  and  particularly  in  figure  8,  there  is  a  strong  correlation. 

From  the  preceding  analysis,  the  following  shape  for  the  dynamic  C.G. 
displacement  can  be  hypothesized:  initially  the  center  of  mass  of  the  dummy 
moves  down  about  1.25"  in  approximately  90  msec,  then  levels  off  or  rebounds 
slightly  as  the  thrust  remains  on.  When  the  thrust  ends,  the  dummy's  center 
of  mass  returns  to  its  initial  position.  This  final  turn  takes  place  at  about 
150  msec,  when  the  displacement  is  roughly  1".  It  appears  that  the  initial 
movement  of  the  C.G.  corresponds  to  the  crushing  of  the  cushion  and  the 
rubber  buttocks  of  the  dummy.  The  displacement  then  levels  off,  or  "bottoms 
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out",  when  the  materials  have  been  completely  crushed.  Finally,  when  the 
thrust  ends,  the  dummy  rebounds  and  the  C.G.  returns  to  its  neutral  position. 
Human  Biodynamic  Simulations 

The  next  step  in  the  development  of  the  model  was  to  analyze  human 
response  data  in  a  manner  similar  to  the  dummy  Tower  Test  data  in  an  attempt 
to  extend  the  model  for  human  occupants.  Because  of  the  paucity  of  human  test 
data  that  would  lend  itself  to  this  type  of  analysis,  and  because  of  the 
extreme  difficulty  in  obtaining  this  type  of  data,  a  different  approach  to  the 
problem  was  deemed  necessary.  Rather  than  relying  on  actual  human  tests  to 
provide  the  necessary  data,  it  was  felt  that  a  biodynamic  gross  motion 
simulation  program  could  provide  the  time  history  of  the  center  of  mass  of  a 
human  occupant  against  which  the  C.G.  Model  could  be  evaluated.  It  must  be 
kept  in  mind  that  when  ejection  seat  tests  are  conducted  with  human  subjects, 
medical  and  physiological  considerations  insure  that  the  tests  are  conducted 
under  ideal  conditions:  the  subjects  are  extremely  well  restrained  and  the 
accelerations  imposed  are  moderate.  The  simulation  program  permits  the 
relaxation  of  these  restrictions  so  that  the  motion  of  the  center  mass  of  the 
occupant  can  be  investigated  for  conditions  with  a  loose  restraint  system, 
higher  acceleration  profiles,  and  different  occupant  initial  positions.  The 
computer  program  used  to  provide  the  dynamic  response  of  various  occupanta 
subjected  to  a  number  of  different  acceleration  profiles  was  the  Calspan 
Simulator  (4),  which  has  been  the  subject  of  several  validation  efforts 
(5,6,7).  The  formulation  of  the  human  occupant  model  has  been  previously 
described  (5).  Briefly,  the  occupant  was  modeled  via  16  segments  and  15 
joints  (figure  21);  the  segment  masses  and  Inertial  properties,  as  well  as  the 
joint  locations,  were  estimated  from  various  sources  (8,9,10).  The  Calspan 
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Simulator  generates  the  complete  linear  and  angular  time  history  of  the  center 
of  mass  of  each  segment  used  and  that  information,  along  with  the  segment 
masses,  allows  the  easy  calculation  of  the  center  of  mass  of  the  entire 
occupant*  By  subtracting  the  instantaneous  location  of  the  occupant's 
C.G.  from  its  initial  position,  the  C.G.  displacement  curve  is  then  easily 
obtained.  As  was  the  case  with  the  Tower  Test  data,  there  are  also  some 
problems  with  this  type  of  analysis.  First  of  all,  the  data  in  question  is 
not  human  test  data  but  only  human  simulation  data  and,  because  of  that,  it  is 
limited  by  the  restrictions  imposed  by  the  simulation  model  that  generated  it. 
The  16  segments  representation  of  the  occupant  may  be  inadequate  to  accurately 
represent  the  C.G.  of  the  entire  occupant  because  the  C.G.  locations  of  the 
individual  segments  were  mere  estimates,  and  errors  in  their  locations  will 
alter  the  position  of  the  occupant's  center  of  mass.  In  addition,  the 
simulation  program  does  not  allow  for  segment  deformation,  nor  does  it  take 
into  account  the  displacement  of  internal  organs.  Finally,  assumptions  had  to 
be  made  about  the  elongation  characteristics  of  the  restraint  system  and  about 
the  deformation  properties  of  the  seat  cushion.  In  spite  of  all  these 
assumptions  and  restrictions,  useful  information  could  still  be  extracted  when 
the  results  from  a  number  of  simulation  runs  were  analyzed  together. 

RESULTS 

Model  Definition 

Having  defined  the  general  shape  of  the  C.G.  displacement  curve,  the 
next  step  was  to  derive  a  model  that  would  exhibit  the  same  behavior.  A 
simple  lumped-mass,  spring  damper  model  was  chosen  to  represent  the  occupant 
and  its  relationship  to  the  seat  (figure  10).  It  was  felt  that  this  simple 
spring-damper  model  could  cover  the  entire  motion  of  the  occupant's  center  of 
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mass  if  the  spring  constant  "k"  (see  equations  in  figure  10)  was  made 
non-linear,  as  shown  in  figure  11.  The  first  slope,  kl ,  represents  the  seat 
cushion  and  the  dummy's  rubber  buttccks  and  remains  in  effect  until  the  C.G. 
bottoms  out.  The  second  slope,  k2,  represents  the  stiffness  of  the  seat  pan 
and  the  seat  structure;  this  constant  comes  into  play  after  the  C.G. 
displacement  has  bottomed  out.  For  simplicity,  the  same  damping  coefficieut 
was  used  for  both  phases. 

Tower  Test  Data  Analysis 

To  test  out  this  model,  9  of  the  best  tower  tests  were  selected  for 
parameter  fitting.  Figures  12  through  20  show  the  seat  acceleration  input, 
the  measured  C.G.  displacement,  the  best  fit  with  a  linear  spring-damper 
model,  and  the  best  fit  with  a  two  piece  non-linear  model  for  these  9  tests. 
The  figures  indicate  that  the  two  piece  non-linear  spring  constant  model  is 
more  accurate  than  the  single  slope,  linear  constant  model,  and  the  fit 
obtained  with  the  non-linear  model  is  quite  accurate.  The  parameters  used  to 
obtain  these  curves  are  shown  in  Table  2.  It  should  be  noted  that  the  heavier 
dummy  requires  a  smaller  bottoming  distance  for  an  accurate  fit.  This  is 
reasonable  since  the  heavier  dummy  will  crush  the  seat  cushion  and  rubber 
bottocks  to  a  greater  degree  prior  to  ejection. 

Human  Biodynamic  Simulation  Data  Analysis 

The  Calspan  Simulator  was  exercised  numerous  times  to  simulate  the 
response  of  human  occupants  of  different  sizes  when  ejected  from  several  of 
the  presently  operational  Navy  aircraft.  The  results  from  four  representative 
runs,  along  with  the  best  linear  and  non-linear  fit,  are  shown  in  figures  22 
through  25.  The  parameters  used  to  obtain  these  curves  are  shown  in  Table  2. 

As  is  evident  from  the  shape  of  the  C.G.  displacement  curves  shown  in 
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figures  22  through  25,  the  behavior  of  the  center  of  mass  of  the  simulated 
human  occupant  differs  from  that  of  the  dummies  tested  on  the  NADC  Tower.  The 
displacement  curves  do  not  exhibit  the  characteristic  trapezoidal  shape  of  the 
dummy  C.C.  displacement,  but  seem  to  be  much  more  linear.  This  observation 
is  further  substantiated  by  the  close  agreement  between  the  linear  and 
non-linear  fit.  The  values  of  the  parameters  used  to  generate  the  curves 
(Table  2)  also  show  that,  whereas  the  dummies  required  a  coefficient  of  .3  of 
critical  damping,  the  simulated  human  occupant  required  a  coefficient  of  .95 
of  critical  damping.  This  is  reasonable  since  the  dummies,  being  composed 
largely  of  metal,  are  much  "stlffer"  than  humans  and  consequently  should 
exhibit  a  less  damped  behavior. 

DISCUSSION 

The  model  presented  clearly  shows  that  it  is  possible  to  obtain  a  close 
approximation  for  the  displacement  of  the  C.G.  of  an  ejection  seat  occupant 
during  a  typical  ejection.  Because  of  limitation  in  the  available  data  the 
model  presented  is  only  applicable  in  the  vertical  or  Z  direction,  but  the 
methodology  outlined  in  this  paper  makes  it  possible  to  evaluate  the  model  in 
other  directions  of  Interest,  and  in  particular  the  forward  or  X  direction. 
Though  one  of  the  primary  reasons  for  undertaking  this  investigation  was  to 
develop  an  analytical  tool  to  improve  the  trajectory  simulation  capability  of 
the  Navy,  the  model  should  prove  valuable  in  the  design  of  rocket  systems  for 
ejection  seats. 
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TABLE  1  -  Selected  Ejection  Tower  Tests  and  Significant  Paraaeters 
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Linear  Model  Non-Llneer  Model 


Occupant 

Weight 

K 

C 

K1 

K2 

Dbot 

C 

(Lbs) 

(Lbs/in) 

(Lbs/ln) 

(Lbs/in) 

(in) 

Tower  Tests 

5841 

145 

1500 

.30 

1000 

2000 

1.25 

.30 

5842 

145 

1500 

.30 

1000 

2000 

1.25 

.30 

5843 

145 

1500 

.30 

1000 

2000 

1.25 

.30 

5844 

145 

150^ 

.30 

1000 

2000 

1,25 

.30 

5845 

212 

1500 

.30 

1000 

2250 

0.75 

.30 

5846 

212 

1500 

.30 

1000 

2250 

0.75 

.30 

5847 

212 

1500 

.30 

1000 

2250 

0.75 

.30 

5849 

212 

1500 

.30 

1000 

2250 

0.75 

.30 

5850 

212 

1500 

.30 

1000 

2250 

0.75 

.30 

Simulation 

Runs 

F18ECQK 

175 

1000 

.95 

1000 

1150 

1.25 

.95 

AV8BEWL 

225 

1150 

,75 

1000 

1250 

1.00 

.95 

F14KDJF 

225 

1000 

.95 

750 

1000 

1.00 

.95 

A4KPAQD 

225 

1000 

.95 

750 

1000 

0.75 

.95 

TABLE  2  -  Results  of  C.G.  Model  Parameter  Fit 
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FIGURE  1  -  Double  Integration  of  Relative  Acceleration  and 
Typical  C.G.  Displacement  Curve 
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FIGURE  4  -  Time  of  the  Two  Turns  for  all  52  Ejection  Tower  Tests 


FIGURE  6  -  Displacement  at  the  1st  Turn  vs  Peak  Seat  Acceleration  for  all  52  Ejection  Tower  Tests 
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FIGURE  7  -  Displacement  at  the  1st  Turn  vs  Maximum  Seat  Velocity  for  all  52  Ejection  Tower  Tests 
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DYNAMIC  C.G.  MODEL 


Ax(f> 

EQUATION  or  motion: 


**■&*+ &*'*(» 

or  6  +  2Y& +  m  Aft) 

where  y  m  JL  u£  *  i 

2m  •  m 

SOLUTION  GENERATES  J  CASES! 

fa)  y  <  <"•  (b)  ?'«»•  (e)  Y  >  “4 

or  Y  =  c«J0  wAere  c</  for  (a),  editor  (b)  and  of  for  ts) 
therefore  6  *  2 cu>0 6  *  u>*6  *  Aft) 

FIGURE  10  -  Dynamic  C.G.  Model  and  Equations  of  Motion 
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FIGURE  11  -  Two  Piece  Non-Linear  Spring  Constant 


-27- 


INVESTIGATION 


I 


28 
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FIGURE  12  -  Monitored  Seat  Acceleration  and  Resulting  C.G.  Displacement  for  Tower  Test  No.  5841 
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Seat  Acceleration  and  resulting  C.G.  Displacement  for  Tower  Test  No.  5847 
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FIGURE  25  -  Input  Seat  Acceleration  and  Resulting  C.G.  Displacement  for  Simulation  Run  A4KPAQD 


